Geometrical characterization of non-Markovianity 
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We introduce a new tool for the quantitative characterisation of the departure form Markovianity 
of a given dynamical process. Our tool can be applied to a generic iV- level system and extended 
straightforwardly to Gaussian continuous- variable systems. It is linked to the change of the volume 
of physical states that are dynamically accessible to a system and provides qualitative expectations 
in agreement with some of the analogous tools proposed so far. We illustrate its prediticve power 
by tackling a few canonical examples. 



PACS numbers: 

The interaction with an environment leads a quantum 
system to dissipate energy and lose its coherence. The 
process, however, needs not be monotonic and the system 
may temporarily recover some of the lost energy and/or 
information. This is the essence of a non-Markovian 
behavior, which can be characterized and quantified in 
many different ways [TJ [2] . One possibility is to look for 
temporary increases of the entanglement shared by the 
system with an isolated ancilla, which amounts to mea- 
sure the deviation from divisibility of the dynamical map 
describing the system's reduced evolution [3] (RHP mea- 
sure). A different approach [4j|5] relies on measuring the 
distinguishability of two optimal initial states that have 
evolved through the same quantum channel, looking for 
any non-monotonicity (BLP measure). Further proposals 
have been put forward, based on the decay rates entering 
the master equation [6], on the Fisher information flow 
[7] , on the use of the quantum mutual information [8] or 
of channel capacities [9j, and on spectral considerations 
[TO] . This variety of tools highlight the multi-faceted 
nature of the problem embodied by a characterisation 
of non-Markovian dynamics and its inherent difficulty, 
which prevents the formulation of a unique tool. 

In this paper, we contribute to the quest for sharp 
tools able to capture the various aspects with which non- 
Markovianity manifests itself and propose a method that 
qualifies non-Markovian evolutions based on the rate of 
change of the volume of accessible states of the evolved 
system. As for the divisibility-based approach [3 , this 
is a characteristic of the map itself that does not de- 
pend on the initial state(s) of the system (nor needs 
to be optimized over them). A quantum evolution is 
Markovian if it is an element of any one-parameter con- 
tinuous completely positive semigroup: in this case the 
process is unidirectional and there is no recovery of en- 
ergy/information/coherence by the system. This im- 
plies that the domain's volume of the dynamical map 



decreases monotonically. On the contrary, we associate 
non-Markovianity of the dynamics to a growth of this 
domain's volume. We thus define a quantifier of non- 
Markovianity as the sum of the (temporary) volume in- 
creases which occur during the time evolution. 

In the case of a single qubit, this can be linked to the 
BLP measure [4 , as the trace distance coincides with 
the Euclidean distance on the Bloch sphere and the pair 
of states that maximize the measure lie on the bound- 
ary of the convex subspace of physical states [TT]. As a 
result, if the trace distance decreases monotonically, so 
does the volume, which is however much easier to evalu- 
ate through the determinant of the dynamical map, as we 
address in details in this paper. Our aim here is to for- 
malize the intuition at the basis of our proposal towards 
the construction of a new tool for the quantitative char- 
acterisation of non-Markovianity. It should be stressed 
that ours is not yet another attempt at the quantification 
of the degree of non-Markovianity of a given dynamical 
process, but the proposal for a novel way to reveal effects 
of an evolution departing from the features of Markovian- 
ity that went so far overlooked. Our proposal enjoys fea- 
tures of practicality and intuition of interpretation that 
are somehow missing from analogous, otherwise equally 
valid quantifiers. 

Systems with finite dimensional Hilbert spaces. Irre- 
spective of the initial open system state, a reduced time 
evolution derived from the unitary dynamics of a larger 
system can always be described by a linear, Hermitian 
map [12], not necessarily completely positive due, e.g., 
to the presence of initial system-environment correlations 
[T3] . A Markovian or memory-less behavior leads to mas- 
ter equations in the Lindblad form p~4j [15] , with the map 
obeying the semigroup composition law. We consider a 
positive trace preserving map 

4> t : p(0) ^ P(t) = MKO)} (1) 
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for the quantum state of a TV-level open system, which 
can be expressed through a generalized Bloch vector 
r, whose components are the expectation values of 
the traceless, hermitian generators of SU(iV), Gi(i = 
1,...,7V 2 - 1), for which Tr[GiGj] = By including 
the identity Go = I/y/N, any state can be written as 
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with r = (l/vJV,r). A systematic construction of 
the {G a } is given in Refs. [H3 [H] and leads us to 

{Ga}a=l l = {Ujk,V jk ,Wi}/V2 with 

Ujk = \j) (M + \k) 01 , v jk = (k\ - \k) 01), 
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where the span of the indices is such that 1 < j < k < 
N, 1 < / < N — 1 and {|m)}^ =1 is an orthonormal basis 
of the open system's Hilbert space. This gives Pauli spin 
operators for TV = 2 and Gell-Mann operators for N = 3. 
Writing the map in Eq. in this basis, one gets 

r t = F(t) f , with F aj3 (t) = Tr[G a </>t[Gp]]. (4) 

As Fop(t) = <5o/3, this is an affine transformation for the 
Bloch vector. Letting qp (t) = Fpo, we have 



F(t) 
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-> r f = A t r + qt/VN. (5) 



The real matrix At can be decomposed as At = 
0\Ti t O\ , where O t "'s are orthogonal matrices and D 
is a positive semi-definite diagonal one. In what follows, 
we will indicate with \M\ the determinant of a matrix 
M. The findings above imply that |F t | = |A t | = |D t |. 
The action of F is given by a first rotation (possibly com- 
posed with an inversion), then a shrink of the Bloch vec- 
tor followed by a final rotation plus a translation. Its 
determinant gives the contraction factor for the volume 
of accessible states, given by the measure of the set of 
evolved Bloch vectors, with respect to its value at t = 0. 

The set of physical Bloch- vectors for an iV-level system 
is given by [18] 
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where aj (r) are the coefficients of the characteristic poly- 
nomial det(xI N — p) with p = j^I^ 
spherical coordinates, the volume element of Bn is 



Ez^i 1 r i G i . In 
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and, by direct substitution, it is straightforward to check 
that any positive trace-preserving map described by 
Eq. Q induces the change 

d N V{t) = \\A t \\d N V{Q). (7) 

Therefore, 1 1 1 1 describes the change in volume of the set 
of states accessible through the evolution of the reduced 
state. In particular, ||A t || decreases monotonically for 
any positive, linear and trace preserving map [15] , and 
so it does for any element of a completely positive con- 



tinuous one-parameter semigroup. Indeed, if (j> t 
with 

£ p = i[p, H]+J2 7«,/3 ( C a pC\ - \{C\C a , p] 



(8) 



7 > 0, and H the Hamiltonian of the system, we get 
| A t | = e _ArTr ^, which is a a constant. 

Time dependent generators of the form in Eq. ([8| with 
j(t) > lead to time-dependent Markovian processes. 
The dynamical map 0t+r,t = exp[T J^ +r £dt], although 
not being part of a dynamical semigroup, is divisible and 
can be written as the composition of two CPT maps 



</>£+r,0 = </>£+r,t 0t,o (Vt,£ > 0). 



(9) 



As a consequence, in this case too the determinant is 
monotonically decreasing. These considerations lead us 
to define a new way to quantify the non- Markovian char- 
acter of a quantum evolution through the variation of the 
volume of accessible states 

dV(t) f d\\F t \ 
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The intuitive meaning of this definition is illustrated in 
Fig. [T] where the time evolution of the determinant is 
explicitly shown for the case of a two-level atom sponta- 
neously decaying in a structured environment (cf. Exam- 
ple 1 for details). The monotonous volume decay char- 
acterizing a Markovian time evolution is contrasted with 
a non-Markovian behavior in which a non-zero accessible 
volume re-appear after being fully decayed. 

Besides the geometric interpretation, such a measure 
has a simple physical meaning based on the change of 
the classical information encoded in the states. Suppose 
that a set of quantum states is given, whose elements 
are characterized by an arbitrary distribution of Bloch 
vectors within B^-, described by a probability density 
p(r). The corresponding differential entropy is 



h\p(v)} 



P( r ) log 2 p(r)dV N . 



(11) 



If such states are taken as the initializations of the map 
4>t, after a time t the probability density function is re- 
scaled as p'(r t ) = p(r t )/\\A\\ and the entropy changes 
accordingly as 



h\p'{v t )]-h\p{v )]=log 2 \\K t 



(12) 
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Thus, a contraction of volume is equivalent to a loss of 
classical information. 

The BLP measure also enjoys an information theoret- 
ical interpretation, but a comparison between the two 
quantifiers is difficult. Indeed, although it is known that 
the optimal states that enter the measure proposed in 
Ref. [4 lie on the boundary of the volume accessible 
throughout the dynamics [11], this does not necessar- 
ily imply a connection with the measure of such volume. 
On the other hand, the difference between Mv an d the 
RHP measure is simpler to describe: as the determinant 
is contractive under composition of positive maps, it fol- 
lows that it does not increase whenever the intermediate 
map 4>t+r,t m Eq. ([9| is non positive. The entanglement- 
based measure of Ref. [3 , on the other hand, is non-zero 
in the less restrictive condition of the map being not com- 
pletely positive. Therefore, if (j)t+ T ,t is positive but not 
completely positive, the map is non-divisible. Nonethe- 
less, we find Mv = 0. 

From a practical viewpoint, the experimental evalu- 
ation of Mv passes through the determination of the 
volume of the set of evolved states found from the con- 
traction of Bn. This is an ellipsoid for the case of a 
qubit. As the number of evolved states that any real- 
istic experimental implementation can sample is finite, 
it would be a precious piece of information to know 
which are the best initial states to use in order to de- 
termine the set of physically accessible states at a given 
time t of the dynamics and its volume. For this pur- 
pose, let us consider N initial Bloch vectors, evolved 
up to time t and arranged as the columns of a matrix 
P t We have from Eq. (J5|, that such vector evolves as 
P t = A t Po + Qt, where the columns of are given by 
the qt's (which provides the form of the evolved state 
for a maximally mixed initial condition). From this it 
follows that |(P t - B t )(P t - B t ) T | = (|A t |) 2 (|P |) 2 . In 
turn, if we choose as initial Bloch vectors the elements of 
any orthogonal basis in M N _1 plus the null vector corre- 
sponding to the maximally mixed state, then their time 
evolutions (arranged to form the matrix P t — Q t ) gives 
the determinant of the map, from which the measure Mv 
easily follows. Therefore, the geometric measure of non- 



Markovianity in Eq. ( 10 ) can be revealed experimentally 
by performing a state tomography at different times for 
N 2 — 1 initial orthogonal states. This will be sufficient 
to evaluate the change in volume of the accessible states 
without the need for prior knowledge about the environ- 
ment or the coupling. For the case of a qubit, this is 
illustrated in Fig. [5] (a), where the three initial Bloch 
vectors corresponding to the canonical basis of R 3 are 
shown to evolve into the extreme points on the semi- axis 
of the ellipsoid that comprises all the possible accessible 
states of the evolution. 

Systems with infinite dimensional Hilbert space. We can 
extend our idea to the less intuitive context of contin- 
uous variable systems, in which the Hilbert space is in- 



finite dimensional and it is not possible to describe a 
state through a finite number of parameters. However, 
the restriction to Gaussian state and Gaussian-preserving 
processes helps overcoming this issue. 

We consider a system made of n bosonic modes k = 
l,..,n, each described by the annihilation and creation 
operators and a\ [corresponding position and mo- 
mentum operators = 1/V2(dk + a\) and pk = 
—i/y/2(a~k — respectively]. Defining the vector of 
operators R = (<?i,pi, . . . , #nPn) T , the commutation re- 
lations can be written as [Rk,Ri] — iftki where Q^i are 
elements of the symplectic matrix ft = with 
u = ( _°x J). A Gaussian state is completely character- 
ized by its first and second statistical moments, given 
respectively by (R) and the covariance matrix a defined 
as 



a kl — -^{{Rk.Rl}) 



(R k )(Ri). 



(13) 



First moments can be adjusted to be null by local unitary 
operations. It can be shown that any evolution resulting 
from the reduction of a symplectic evolution on a larger 
Hilbert space can be described, in terms of second mo- 
ments, by the equation 



a t -+ Xi a X t + Y t 



(14) 



where X and Y are 2n x 2n real matrices fulfilling the 
relation Y + id — iX T QX > 0. Viceversa, any evolu- 
tion of this kind may be interpreted as the reduction of 
a larger symplectic evolution [20 . In full analogy with 
what we did in Q, by choosing a basis {Gj} for the space 
of 2n x 2n matrices, we can write such map as 

a * = ^[X?GkX t Gj]Tr[aoGk]Gj+Tr[Y t Gj]Gj. (15) 

jk 



Eq. (15) can be recast as the R 4Ar affine transformation 
so —> s t = X t So + Yt. We then define a measure of non- 
Markovianity in a way fully analogous to what has been 
done above for the case of a discrete- variable system, i.e. 



as in Eq. (10) with th replacement F t —> Xt- 



If the evolution is unitary, the associated transforma- 
tion has a constant determinant, equal to one. For a 
single mode in a generic noisy Markovian channel, we 
have EE] 



—Ft 

a t = e a 



(l-e-^Voo, (16) 



which gives |X^| = e _4Ft . Therefore, every increase in 
| X t | signals non-Markovianity. 

Having introduced our formalism, we now illustrate 
our proposal with the aid of a few significant examples. 
Example 1: Spontaneous emission into a leaky cavity. 
Consider a single two-level atom with transition fre- 
quency ujq interacting with a vacuum electromagnetic 



(a) 



(b) 




FIG. 1: (Color online) Time evolution of the determinant 
of the map for a generic Markovian [panel (a)] and non- 
Markovian dynamics [panel (b)]. Quantitatively, the curves 
displayed in the figure correspond to the spontaneous emission 
of a two-level system in a resonant leaky cavity (cf. Exam- 
ple 1) with the Lorentzian spectral density given in Eq. (17). 
Panel (a) shows the (Markovian) case with 70/A = 0.1 (bad 
cavity limit), while panel (b) shows a (non- Markovian) evolu- 
tion in the good cavity limit with 70/A = 10. We also picture 
the set of accessible states, whose volume changes in time 
according to the behavior of the determinant. 



field having a Lorentzian spectral density (mimicking a 
leaky cavity) [5]. Taking 



J(w) 
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2tt (w - A - w) 2 + A 2 ' 
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where A is the detuning between the atomic and the 
cavity frequency, the atomic state at time t reads 



PA(t) 



r(t)lW + 
r(t)Vo + 
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with fl± = A - iX ± a/(A - iA) 2 + 270A and 
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The evolution of the Bloch vector is ruled by 

9ter(t) 3mT(t) 
A f = I -3mT(t) 9ter(t) 





(18) 



(19) 



(20) 



whose determinant, |A^| = |T(t)| 4 , is shown in Fig. [T] 
against the dimensionless time Xt. The corresponding 
non-Markovianity measure My is reported in Fig. |2j from 
which it is clear that a strongly non-Markovian behavior 
is found for a resonant coupling and invoking the so- 
called good-cavity limit 70 ^> A. A similar result is ob- 
tained with the RHP measure, [3] which is given by the 
integral of (l/2)9te[dt In V{t)\. This is in agreement with 
BLP as well, which turns out to depend on \T(t)\ [5]. 
Example 2: Pure dephasing. Let us consider a qubit 
undergoing a purely dephasing dynamics, expressed in 
terms of a decoherence factor v(t) as 
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FIG. 2: (Color online) Panle (a) Example of an ellipsoid rep- 
resenting the image of a random map for a qubit. The canon- 
ical basis of R 3 is mapped onto the corresponding colored cir- 
cles highlighted on the ellipsoid, from which the volume can 
be obtained. Panel (b) shows the non-Markovianity mea- 
sure J\fv as a function of the detuning, for the spontaneous 
emission dynamics into a reservoir with Lorentzian spectral 
density (cf. Example 1). 



In this case, the BLP and RHP measures coincide [22], [23]. 
In fact, the trace distance optimised over the initial qubit 
preparation gives us D\p°^ 1 , p^ 1 ] = \v(t)\. In turn, \v{t)\ 
is exactly the value of the concurrence between a sys- 
tem qubit and ancilla, initially prepared in a maximally 
entangled state and undergoing a unilateral dephasing 
mechanism. As for our proposal, the evolution of the 
Bloch vector is determined by the matrix 



A t = -3mv(t) 9tei/(t) 
V 1. 



(22) 



which has ||At|| = \v(t)\ 2 . The geometric measure My 
thus gives the same behavior predicted by the other two 
quantifiers. 

We have proposed a geometrically motivated quanti- 
fiers of non-Markovianity that is explicitly linked to the 
variations in the volume of the physical states dynami- 
cally accessible by a given open open system. From an 
information theoretical perspective, such measure, which 
provides predictions that are, in general, in qualitative 
agreement with those coming from some of the most pop- 
ular tools for the characterisation of non-Markovianity 
proposed to date, is linked to the loss/regain of classical 
information oven the evolving system. We have shown 
how an estimate of the proposed measure is possible 
through only a polynomial number of measures, while 
the proposed formal quantifier itself enjoys a straight- 
forward extension to the Gaussian continuous-variable 
scenario. We have illustrated our proposal through a 
series of examples. We hope that the appealing aspects 
of practicality and intuitive nature of our proposal will 
soon spur the attention of the community interested in 
open-system dynamics. 
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